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Abstract 


Model for studying coupling dependence on lattice spacing a in gluo- 
dynamics is suggested. The model predicts g ^ go > 0 with o ^ 0. Free 
energy density in the model does not depend on temperature. 


1 Introduction 

As it is known, in lattice glnodynamics only the closed surfaces formed by 


plaquetts 


□x/ii/ = Ufj, (x) Uv [x + y) {x + v) ul (x ); (x) € SU (N) (1) 


contribute to the expansion in powers of (3 of the partition function 



where g is the coupling constant. A surface may be self-intersecting, but it 
should not intersect another, since in this case those two surfaces must be treated 
as one. Since no additional restrictions are imposed, each surface may be treated 
as a random membrane. 

There is some space in and out of each surface which is unavailable for 
another surfaces. We prescribe to any surface Sk some effective volume 14 = 
V {Sk), that depends both on area and shape of Sk. In fact it is a volume 
of a minimal imaginary external shell, that films over each plaquette of the 
considered surface with the layer of effective thickness b no less than 1/2 link. 
For instance, surface of a cube with n links on edge, has the volume of 6 x 
{n + l)^xl/2. One must also take into account some effective redundant volume 
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V (S) — bS, that appears because a gap between the surfaces can’t be reduced 
to one link distance because their shapes are not rectangular in general. 

Since the states, which differ by of surface shifts give the same contribution, 
the number of states contributed by a surface must be proportional to the 
volume available for such shifts. Therefore, the pattern of exactly solvable model 
PP may be adjusted to develop a model intended to compute mentioned closed 
surface contributions in lattice gluodynamics. Since within model framework all 
orders in /3 are taken into account, model application region may be extended 
to the weak coupling area g < 1- 


2 Model 


Let the four dimensional volume Vtot = contain r surfaces with volumes 

Vk = V (Sk), k = 1 ,..., r. The available volume for any surface is Vtot — 

Let 5 R 5 be the number of configurations, taken by an isolated shape S without 
shifts, then total number of configurations may be computed ’in the spirit’ of a 
van der Vaals approximation 


where |2] 


and 


CX> 


^ = EE 

r=0{S,} 


{Vtot 




■ exp 



li^s, 

k 



X > 0 

X <0 


r! 

2TTi 


n-\-0-\-ioC 
’ -1-0 —200 


gP^p-l-r- 


dp 


f In (g-yN) for SU {N) ■, N>2 

I lng-2 for SU{2) 


( 3 ) 

( 4 ) 

( 5 ) 


Further we shall change the order of summation and integration. To make it 
safe we shift the integration path 


{Vtot - XI Vk)l 


1 fC+ioo 

_ / QP{Vt‘>t-T.lVk) p-l-r 

27rz J c—200 


( 6 ) 


in a such a way that c > Repa for singularity p = Pa with the largest Rep and 
for the partition function we may write 


Z = 


27rz 


^C-|-200 


/'C-\-ioo 


cPVto 


aPVto 


exp {tSk - pVk} ^Skdp 


r—0 


{Sj} k 


2TTi 




( 7 ) 


r=0 


where 


3?= ^ exp{tS-pV{S)}^sdS 


( 8 ) 


S=Srr 
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Further we replace summation in by integration 


pCX) 

SR= ^ exp{tS -pV{S)}^sdS ^ / exp{tS -pV{S)}disdS (9) 
S=Sm 

Possible consequences of such approximation are discussed in Appendix I. 
Having collected everything we may finally write 


^c+ioo 


z = 


expjpVtot} 


1 


27ri Jc-ioo P-^ 


dp 


( 10 ) 


With infinitely increasing lattice volume Vtot = partition function 

may be with good accuracy written as 


Z ~ expjpa (t) Vtot} 


( 11 ) 


where Pa is the singularity of integrand in m with the largest Rep. 

There two possible regimes. The first one is realized when the solution of 


p — 5R = 0, 


( 12 ) 


dominates giving the rightmost singularity in a complex p-plane. Another 
regime is realized when singularity of K proper dominates. 

Note, that in fact we must include the self volume of Sj into the available 

r 

one, i.e. replace {Vot - Yll by {Vot - 14 + Vj). A simple but bulky 

1=1 

computation allows us to write for the partition function 

1 exp {pVtot} f oso/o 1 j no\ 

v-’K ( 13 ) 

However, this correction doesn’t change the results appreciably, since the posi¬ 
tion of singularity of d^/dp is the same as that of 5R. 

Domination and possible interchanges of regimes are defined not only by 
specific form of 3?, but are also determined by the choice of the coupling con¬ 
stant (and consequently t) dependence on lattice spacing a. Indeed, free energy 
density is 

T T 

F = -—\nZ -—VtotPa (t) = -a~'^pa (t) (14) 

where V = is a veritable volume of the system and T = {Nra)~^ is the 

temperature. For any pa which does not vanish with a —> 0 as does x const, 

free energy density becomes infinite, so such state becomes unavailable for the 
system. Therefore, we claim 


Pa {t) = a^Pa (t) 


(15) 
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where Pa (t) is assumed to be finite. Note that when dependence g on a is 
specified and a ^ 0, free energy density become constant F = — lima^o Pa (t) = 
—Po- In particular, it doesn’t depend on temperature. 

Another two basic assumptions of the model are 


(Us ^Cexp{pS} 

(16) 

(see Appendix II) and 

V (S) -bS+ (S/af 

(17) 

Some reasons for above assumption^ are given in Appendix III. 
Therefore we may finally write 


pOO 

K = 3?(t,p)~C' / exp {— rS* — p (5'/cr)“} 

-t Sm 

(18) 

with 

T = pb — p — t 

(19) 

Integral in 11811 diverges for all r < 0 and converges for all r > 0. In 
point T = 0 integral in 11811 converges to finite value 

a singular 

3?(0,p) ~ —p~iV (-,p{Sm/<j)°‘) 

a \a / 

(20) 


for any Rep > 0. If singularity r = 0 of K(r,p) is leading in a limit a —> 0, it 
specifies the dependence t on a. Indeed, from pa (t) b = a'^Pa {t)b = p + ta we 
get -where Pq = lima^o-Pa (t), hence > 0 with 

a ^ 0. 

Since d^/dr < 0 for r > 0, function 3?(t,p) is monotonously decreasing 
in this area. Therefore, if the solution pa (t) = "Silpa (t) b — p — t,pa (t)] of (IT^ 
exists, it is unique. There is some range of values of t and a where such solution 
exists and varying those parameters we may move Pa (t) in this area. When 
dependence t on o is specified, i.e. t = ta this area shrinks into line Pa (ta) 
defined by the single parameter a, and according to (Cnj Pa (ta) 0 with 
a ^ 0. 

If ta is chosen in a such way that limQ_,o ta = to = const, then g^ —> gg 
= > 0. The case of ta ^ to ^ —oo is unacceptable, since it means 

g ^ oo with a ^ 0. It is easy to check that we cannot choose g ^ 0 (ta ^ to ^ 
+CX)) with a —> 0 because there is no solution p = 3? (t,p) —> 0 for such choice 
of ta- Indeed, in this case from we obtain 

^{Ta,Pa)^Cy^ ^ ’ / exp{-T„5}5^+“"-ld5 (21) 

or 

(fi{Ta,Pa) ^ Ct--^ £ r(7 + aj,rA) ^_,, (-r-“p)” (22) 

_ r^O 

^As it follows from 0, role of specific form of V (s) with a ^ 0 is essentially diminished 
due condition 1151 . 
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that with 0 


r (A,a;) 


N-l 

n—0 


r(l-A + n) 1 
r(i-A) (-x)- 



(23) 


leads to 

3?(ra,p) ~ -CS2, {p + ta)~^ exp {{fj. + ta) Sm) (24) 

so Pa = ?fi(Ta,Pa) ^ —c)o with 0 ^ 0 , but we claim pa ^ 0 in such a limit. 
Hence, there is no asymptotic freedom in suggested model. 

The possibility for QCD (and gluodynamics in particular) to be a non- 
asymptotically free theory have been discussed for years. There are reasons 
to believe, that QCD is not perturbative at a ~ 0 0. On the basis of today’s 
numeric computations, it is difficult to anticipate the behavior of g in the limit 
of a —> 0, taking perturbative calculations as a guidance. Moreover, numerical 
studies IS] showed deviations of the Callan-Symanzik /3-function (3cs{g) from 
perturbative result when the correlation length begins to grow. These devia¬ 
tions are of such a pattern, as if the theory approaches the fixed point at 
which f3cs = 0 and consequently the theory is not asymptotically free. Solid 
arguments in favor of such behavior of Pcs (go) were given in [S]. Data on 
deep inelastic scattering does not eliminate the fixed point |Zj. Phenomenolog¬ 
ical analysis of available monte-carlo lattice data in the S'C/(2)-gluodynamics 
shows no contradiction with the fixed point of Pcsig) located at go ~ 0.563 
|S]. Analytical estimations |S] also favor go 0. 


3 Conclusions 

The pattern of exactly solvable model P is adjusted to develop the random 
membrane model for lattice gluodynamics. We make use of the fact that only 
closed surfaces formed by plaquetts contribute to the expansion of the lattice 
gluodynamics partition function in P powers. Since within model framework all 
orders in P are taken into account, model application region may be extended 
to the weak coupling area g < 1. Arguments for the main assumptions, m 
and 03 of the model are given in Appendix II and Appendix III. 

Model predicts 5 —> 5 o > 0 with a ^ 0 and independence of free energy 
density on temperature. 


4 Appendix I. Discrete Volume 

For the discrete variable m one may define 


( rrf 

for 

m > 0 


m\={ <55 

for 

m = 0 

(25) 

1 0 

for 

m < 0 
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and write a series 


(26) 


oo oo 



n— — OQ n—0 


which is a discrete version of Laplace transform, called Z - transformation (see 
e.g. [ID]). The inverse transform 



(27) 


WC-ZTT 

may be used instead of 0 , so for partition function we obtain 



(28) 


that leads to 



(29) 


where 


SR = ^ exp {pV {S) + tS} Ks 


(30) 


s 


and one may finally write 



(31) 


Since exp{5R —p} is entire function of SR — p, singularities of the integrand in 
m and in m are located at the same position, and the difference between 
corresponding expressions for partition functions disappears with Vtot —>■ oo. 

5 Appendix II 

For simplicity we consider here the closed two-dimensional surface in the three- 
dimensional space, instead of four-dimensional space in suggested model. The 
slice between planes x = Xc and a; = a;c -f 1 is a two-dimensional surface closed 
in the x-direction. Its borders are located in the above planes. The borders are 
equal and are closed loops of length L^- The area of this slice is x 1. For N 
slices of common area Syz we get 




(32) 
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or 


where 


{Syz) = tt- / expIpR^z}^^' 
Sttz 

_ pOO 

R = Y^ Rlc-p^ ~ / RLe-P^dLdp 

T 'Z Z/min 


Summing over TV we finally get 




iV=l 


It is known (see e.g. CH) that 

Rl = CL^~^ exp {mL} 
where C, A, m are constants ^ and we get 

poo 

R‘ 2 iC exp {—{p — m) L} dL 

" Z/min 


= [p-m) '' cr (A, (p - m) Lmin) cz {p - m) ^ CT (A). 
Having written 


(33) 

(34) 


(35) 


(36) 


(37) 


1 - i? = 1 - CT (A) (p - m)~^ 


^ (n + A)! 

^^(n + l)!r(A) 


(-x)"(Cr(A)) 


(n+l)/A 


(38) 


where x = p — m — (CT (A))^^'^,we come to an integrand in (I35II having a simple 
pole in p = /i = m + (CT (A))^^^. Since A > 0 we get p = p > m, so singular 
point of i? at p = TO is located leftward in complex p-plane. Therefore, the pole 
is a leading singularity and for Syz 1 we may write 


(Syz) — exp {pSyz} (39) 

If one assumes that 

5R (Sx) ^ exp {pSx} (40) 

where Sx = S — Syz part of surface that consists of x-planes plaquetts, one 
comes to 

3? (S') ~ / exp{/i (S — Sjf)} exp{pSx} dSx — Sexp {/iS} (41) 

•/ Smin 

Unfortunately such simple estimation allows to exclude corrections neither in a 
power type factor, nor in p. Nonetheless, even in this case expression for 5ft (S) 
will not contradict to the assumption 

^For closed loops without intersections C ~ 6/5;A ~ 4/3;m ~ 1 1111 . but intersection 
doesn’t change those constants drastically. 
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6 Appendix III. Packing of surfaces. 

The proper effective volume of the surface S is very close to bS^ but with increas¬ 
ing number of surfaces and their areas the packing problem appears. Despite 
four century history, this problem has been more or less solved only for objects 
of simple form (mainly for spherical ones) ^|E|- In particular, for ’random’ 
balls packing in a spherical bag of volume v, one can define the density p = vb/v 
where vb is the volume of the balls. The density depends on many conditions, 
but roughly it may be computed as m 

(42) 

where N is the number of balls. One may expect similar density behavior in a 
case of random volume of balls. Let v = vb/N is an average volume of a ball, 
hence one may write 

'P p 3 1 — 

-n-—T = (43) 

^-^(vb/v) 2 

If for fixed average volume v we increase volume of balls us, in accordance 
with 63 , by so doing we increase v. One may expect, that similar relation 

V^ciVb+C2V^ + .... (44) 

is true for a single ball volume Vb and volume V which this ball effectively 
occupied. Taking this relation as a pattern, we assume that there exist some 
constants b > 0; a > 0 and 0 < a < 1, such that for objects of arbitrary form 

V (S) -bS+ {S/a)°‘ (45) 

Although expression 63 is regarded only as the model assumption, there is 
at least one more argument in favor 63 . For the convex body the volume of 
parallel shell with thickness 6 is {S) = Sb + Mb^ + ^b^ where M is an average 
curvature integral na Dimensional method allows to assume that M ^ 
that corresponds a = 1/2 in 11711 and allows to expect that with varying a in 
63 a reasonable description may be found for a body of arbitrary form at least 
for S' ^ 1. 
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